Short range stationary patterns and long range disorder in an evolution 
equation for one-dimensional interfaces 
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A novel local evolution equation for one-dimensional interfaces is derived in the context of erosion 
by ion beam sputtering. We present numerical simulations of this equation which show interrupted 
coarsening in which an ordered cell pattern develops with constant wavelength and amplitude at 
intermediate distances, while the profile is disordered and rough at larger distances. Moreover, for 
a wide range of parameters the lateral extent of ordered domains ranges up to tens of cells. This 
behavior is new in the context of dynamics of surfaces or interfaces with morphological instabilities. 
We also provide analytical estimates for the stationary pattern wavelength and mean growth velocity. 

PACS numbers: 47.54.-r, 68.35.Ct, 05.45.-a, 79.20.Rf, surfaces 



Pattern formation is ubiquitous in nature and one 
of the most fascinating features of nonequilibrium 
systems [l[ . Typical examples can be found in many 
startlingly similar interfaces which emerge in very 
different processes, such as growth of amorphous [2| 
and epitaxial thin films Q, or erosion by ion beam 
sputtering (IBS) These pattern forming sur- 

faces can be classified into different categories ac- 
cording to the stationary or time-dependent (coars- 
ening) behavior of the typical pattern length scale 
I. Actually, a large effort has been devoted recently 
to describe these systems through height equations 
since these provide compact and efficient ana- 
lytical/numerical descriptions that successfully de- 
scribe global morphological properties such as ki- 
netic roughening, and surface pattern formation and 
coarsening. Specifically, in order to assess the pre- 
dictive power of (Id) height equations, it would be 
important to produce criteria for the presence or ab- 
sence of coarsening. In Ref. [|| up to four different 
scenarios have been proposed depending on the be- 
havior of / as a function of the pattern amplitude A. 
The conclusion is that coarsening stops (interrupted 
coarsening) if the function 1(A) attains a maximum, 
after which the amplitude increases indefinitely with 
time. A similar conclusion was reached at by Krug in 
Ref. Q, where it is actually conjectured that no (Id) 
local height equation can "describe the emergence 
and evolution of patterns with constant wavelength 
and amplitude" . Hence, a Id counterexample of an 
interface equation [j| leading to a stationary pat- 
tern with constant wavelength and amplitude would 
be interesting to improve our understanding of in- 
terrupted coarsening and, specifically, of the type of 
nonlinearities that induce it Q . 

Recently, a continuum 2d model has been intro- 
duced which describes interesting (sub)micrometric 
features of surfaces eroded by IBS hj. The re- 



sulting interface equation was studied numerically 
for moderate system sizes suggesting the occurrence 
of a stationary ordered pattern and interrupted 
coarsening [ll|. Here, we derive the Id counterpart 
of the height equation in [Icl ] from a physical model 
of IBS and perform a systematic numerical analysis 
of its coarsening properties in relation with Krug's 
conjecture. 

Following the standard assumption made in hy- 
drodynamic models of acolian sand dunes we 
will consider that ripples formed under IBS are 
translationally invariant in the y direction; addition- 
ally, we assume symmetry under x — > —x, as occurs 
under normal incidence conditions for the bombard- 
ing ions [l3|. As proposed in [l(|, the evolution of 
the thickness of the mobile surface adatoms layer R 
and the height of the bombarded surface h is pro- 
vided by a pair of coupled equations, namely, 
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where T ex and T a( i arc, respectively, rates of atom 
excavation from and addition to the immobile bulk, 
(1 — (ft) measures the fraction of eroded atoms that 
become mobile, and the third term in Eq. |T]) de- 
scribes thermal diffusion of mobile adatoms. 

The rate at which material is sputtered from the 
bulk is described by microscopic derivations 14| and 
depends on the local morphology of the surface 



T e , = a [l + a 2 dlh + a 3 (d x h) 2 ]. 



(3) 



where ao is the sputtering rate for a planar surface. 
The rate of nucleation is also related to the local 
shape of the surface, and is given by 

Tad = 70 [R(l + I2&lh) - Reg] , (4) 

where R eq is the thickness due to the mobile atoms 
that are thermally generated even in the absence of 
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bombardment, and 7^ is the average time between 
nucleation events. After a multiple scale expansion 
of H])-([l]), ^ can be adiabatically eliminated from 
Eqs. (HJ) and ([5]) [HQlj, obtaining, to lowest order 
near threshold of the morphological instability, the 
following Id equation for the evolution of the surface 
height, 

8th(x,t) = -isdlh-Kdih + \ 1 (d x h) 2 -\ 2 dl(d x h) 2 1 

(5) 

where parameters are related to those in (HJ)— (J4J) and 
depend on the experimental conditions. To the best 
of our knowledge, the deterministic Eq. (JSJ) has not 
been systematically studied. We will restrict our- 
selves to positive values of v and K, which are re- 
quired in order to produce a long-wavelength insta- 
bility. Moreover, Ai and A2 are required to have the 
same sign for mathematical well-posedncss, as shown 
in [l5], [l|| ■ If the signs of the non- linear terms are 
simultaneously changed, Eq. (|5|) remains invariant 
after h — > —h. Thus, we will only consider positive 
values of these parameters. For A2 = 0, Eq. ([5j re- 
duces to the celebrated Kuramoto-Sivashinsky (KS) 
equation (l7l . [l8j . which is a paradigm of spatio- 
temporal chaos. Its nonlinear term stabilizes the 
system and a (disordered) pattern develops that is 
characterized by a wavelength that does not coarsen, 
and by chaotic cell dynamics. On large length scales, 
the KS system can be effectively described by the 
stochastic Kardar-Parisi-Zhang (KPZ) equation 
paradigmatic of kinetic roughening. In particular, 
the surface roughness (global rms width, W) [2(| for 
a KPZ interface scales as a power law with the lat- 
eral system size L. On the other hand, for A2 7^ 
and Ai = 0, Eq. (|5|) reduces to the "conserved" KS 
equation. This equation has been studied in the con- 
text of amorphous thin film growth [l5T ] and step 
dynamics on vicinal surfaces [211 ]; in this case, the 
linear instability evolves into an ordered pattern of 
paraboloids with uninterrupted coarsening. 

In order to reduce the number of parameters and 
simplify the analysis of ((5]), we rescale x, t and h 
by (K/v) 1 / 2 , K/v 2 and v/\\, respectively, resulting 
into a single-parameter equation, namely, 

d t h{t,x) = -dlh-d*h+{d x hf-rdl{d x h) 2 , (6) 

where r = (v\i)/{K\x) is the (squared) ratio of a 
linear crossover lengthscale to a non-linear crossover 
lengthscale. We have performed a numerical in- 
tegration of © using a fourth-order Runge-Kutta 
method and the improved spatial discretization in- 
troduced by Lam and Shin [22j for the nonlinear 
terms. We have used periodic boundary conditions, 
lattice constant Ax = 0.5 and time step At = 0.01, 
checking that results do not differ significantly for 
smaller space and time steps. The initial height val- 
ues were chosen uniformly distributed between and 
1 and statistical data were obtained as averages over 
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FIG. 1: Height profiles for Eq. @ with r = 50 on a 
system size L — 128 at times (a) t — 2, 6, 11; (b) t = 11, 
24, 41; (c) t = 41, 70, 103; (d) t = 103, 153, 205; (e) t = 
205, 302, 425; (f) t = 425, 764, 1024. Height profiles at 
different times evolve by increasing the maximum value 
of h with time. All units are arbitrary. 



250 random initial conditions. The standard system 
size of our simulations was L = 512 (1024 nodes), 
except when other values are indicated, and the pa- 
rameters fixed to v = 1, K = 1, Ai = 0.1, varying 
A2 in order to check for the different values of r. 

In Fig. [1] the evolution of the height profile is 
depicted for r = 50. Starting from an initial ran- 
dom distribution, a periodic surface structure with 
a wavelength of about the maximum of the linear 
dispersion relation [l[, namely kmear = 2y/2n, arises 
and the amplitude of h increases (Fig. QJi) . At later 
stages, the conserved KPZ nonlincarity d 2 (d x h) 2 
(cKPZ) induces coarsening of the ordered cell-like 
structure, wherein the cells grow in width and height 
and the number of cells decreases. This coarsening 
is such that smaller cells are "eaten" by larger neigh- 
bors (Figs.QJi-e). We show in Fig.[5^i the time evolu- 
tion of the mean height hi,(t) = 1/L^ a , h(x,i), the 
wavelength l(t), and the amplitude A(t) of the pat- 
tern defined as the mean lateral distance between 
two consecutive local minima and the mean verti- 
cal distance from a local minimum to the next lo- 
cal maximum, respectively. As seen in Fig. [2ja, for 
t > 70, the non-conservedKPZ term (d x h) 2 becomes 
relevant and the mean height of the surface hr, starts 
increasing to reach a constant velocity (Fig. [If). At 
the same time, the coarsening process slows down 
until stopping completely in the stationary state. 
This behavior suggests that the cKPZ term, which 
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(a) (b) 

FIG. 2: (a) Time evolution of mean surface height av- 
erage h,L(t), wavelength l(t), amplitude A(t), and global 
width W(t), for r — 50; (b) corresponding surface struc- 
ture factor as a function of wave number q at times 
t = 4, 15, 30, 60, 125, 250, 1129, and 10653, bottom to top 
(curves are offset vertically). All units are arbitrary. 




( a ) (b) 

FIG. 3: (a) Height profile at t = 15000 for a system size 
L = 10000 and r = 50; (b) local width vs window size 
at t = 15000 for several values of r and L = 8192. Inset: 
same data plotted vs xo/(2V&r). All units are arbitrary. 



acts at small scales, induces the order and the coars- 
ening process until the slopes and the characteristic 
wavelength of the pattern are large enough to make 
the KPZ term no longer negligible. At this time, 
the KPZ term interrupts the coarsening process, as 
claimed in [1, [l(| , and a constant average velocity 
value is achieved as a consequence of a constant av- 
erage of slopes across the interface. Thus, the fi- 
nal wavelength of the pattern depends on the inter- 
play between the cKPZ and the KPZ terms. Fig. [2k 
also shows the global surface rms width or roughness 
W(t) [2(|, and Fig. [2b shows the surface structure 

factor, defined as S(q,t) = (ji(q,t)h(—q,t)) , where 

h(q,t) is the Fourier transform of the field h(x,t). 

The first peak in the structure factor indicates the 
dominant wavelength of the pattern. We can see how 
this peak moves to larger wavelengths (coarsening) 
until a fixed mode is reached which corresponds to 
the stationary value I « 36 in Fig. [2k- At this time 
coarsening interrupts and the amplitude saturates to 
a constant value. Nevertheless, for large enough L 
and at long distances, the profile disorders in heights 
(Fig. [3k) , although the lateral cell- like order is still 
preserved for intermediate distances (Fig. [3k, inset; 
note the difference in scales between the x and y 
axes). This disorder reflects in the power law be- 
havior of S(q,t) for q much smaller than 2ir/l and 
long enough times (Fig. [2Jd) or, equiyalently, in the 
behavior of the local width w(xq) [2fJ] at long times, 
displayed in Fig. [3b function of window size, 
xq. Due to the parabolic shape of the cells, the local 
width scales as w ~ Xq for distances smaller than 
the cell size, reaching a plateau for intermediate dis- 
tances, and finally increasing as a (smaller) power of 
Xq for large enough distances (kinetic roughening). 
The plateau in w(xq) is related to the lateral order of 
the pattern, reaching up to several tens of cells (see 
e.g. in Fig. 0k an ordered domain containing over 30 
cells for r = 10). The effective exponent character- 
izing the long distance behavior of w(xq) increases 



towards its KS value (1/2) for decreasing r. Indeed, 
the profile is more disordered for small r, as can be 
observed from Fig. 0k, to the extent that for r w 0.2, 
and independently of L, the secondary peaks in the 
structure factor vanish completely (not shown), and 
only a weak peak about the linear instability per- 
sists, as in the KS equation. Thus, the KPZ term is 
seen to act at larger scales and is responsible for the 
disorder of the profile, while the cKPZ terms domi- 
nates at smaller scales with a trend to order the cells 
vertically. 

As indicated above, for large r values and after the 
initial times, the cell- like structure is well ordered 
and the length scales of the profile are large enough 
so that the linear fourth order derivative can be ef- 
fectively neglected. Thus, we can rescale x — > r 1 / 2 ^, 
t — > rt and h — ► h to obtain an effective parameter- 
free equation. This means that, for large values of 
r, the solution of Eq. ([6]) remains unchanged if we 
rescale lengths by r 1 / 2 and times by r. In order to 
check this hypothesis we present in Fig. 0Jd the fi- 
nal amplitude and wavelength of the structure for 
different values of r. We can observe that, for large 
values of r, the height amplitude does not change 
with r and the wavelength I is proportional to r 1 / 2 
as predicted above. This behavior can be also seen 
in the inset of Fig. (3}dj where we rescale the horizon- 
tal axis by I = 2v6r (see below), obtaining almost 
perfect collapse of the curves for large r. 

Further progress can be made by recalling that 
the solution of ||SJ) with Ai = is a periodic juxta- 
position of parabolas of the form [l5|, [2l| 

4A 

h(x)=A- 1F x 2 . (7) 

We can consider this function as an approximate so- 
lution of Eq. ([6|) for large values of r. In Fig. [5k 
we compare the numerical profile with the function 
(0 for r = 50. Using this solution and assuming 
as a condition for the final (interrupted) structure 
that the averaged nonlinear contributions along one 
period must be equal, we obtain a relation between 
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FIG. 4: (a) Height profiles at t = 15000 for (top to 
bottom) r = 0.1, 0.5, 1, 5, 10, 50, and 100; (b) stationary 
wavelength I and amplitude A as functions of r (statisti- 
cal errors are smaller than the symbol sizes) . The dashed 
line corresponds to I = 2V6r. All units are arbitrary. 




FIG. 5: (a) Stationary height profile at t = 15000 for r = 
50. The solid line represents the approximate solution 
given by (J7J) for I — 2s/&r and A — 7; (b) mean height 
for different values of r as a function of t/r. The dashed 
line is given by p|. All units are arbitrary. 



the final wavelength and r, namely I — 2y6r. This 
function is represented in Fig. 0}d, fitting accurately 
the stationary wavelengths obtained numerically for 
large values of r. The net mean growth velocity of 
the height average is only due to the KPZ nonlinear 

term and is given by v 



Xl[d x h(a 



Assuming 



(0) as an approximate solution and integrating over 
one period, we obtain the net mean profile evolution. 
It reads 



r-independent if we rescale t — > rt, as already indi- 
cated. Furthermore, at long times, the growth veloc- 
ity given by ([H]) agrees accurately with the numerical 
observations for a wide range in r. 

In conclusion, we have derived a novel determin- 
istic Id equation with shift symmetry in the context 
of ion-beam sputtering and performed a numerical 
analysis to obtain some relevant information about 
its solutions. For large values of r we have estimated 
the final pattern wavelength assuming a parabolic 
solution analogous to that of the "conserved" KS 
equation and checked our assumption by compar- 
ing with the numerical mean height evolution. The 
resulting single-parameter Eq. ([6]) interpolates be- 
tween the KS equation for r = 0, which presents a 
chaotic solution and no coarsening, and the "con- 
served" KS equation for r — > oo (upon rcscaling 
h — ► h/r), which displays unbounded coarsening. 
This behavior is similar to the Id convective Cahn- 
Hilliard (cCH) equation studied in (23|. However, 
as reported in Ref. (23[, coarsening docs not in- 
terrupt in the cCH system for a whole range of 
parameter values but, rather, proceeds logarithmi- 
cally for long times, whereas Eq. ([6]) does present 
interrupted coarsening with an ordered pattern of 
constant wavelength and amplitude for intermediate 
distances. This behavior provides a new scenario for 
the classification of surface coarsening phenomena 
in Ref. [|| , where all the studied evolution equations 
display perpetual coarsening, or else develop a pat- 
tern with a frozen wavelength while the amplitude 
continues growing without bound in the course of 
time. Regarding Krug's conjecture Q, the present 
example reinforces its validity for long range order 
properties, since the pattern produced by Eq. 
disorders at large length scales. However, our results 
suggest that it is still possible to stabilize a well or- 
dered pattern over (intermediate) distances ranging 
even up to several tens of cell sizes in contrast with 
other equations in this context. 
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Ai— t, 
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where we have substituted I = 2v / 6r, and A f» 7 
is obtained from Fig. [4Jd. As seen from Fig. 03, 
for large values of r the evolution of Jil becomes 
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